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ABSTRACT 

This  document  describes  the  design  and  Implementation  of  the  "fly-by- 
wire"  despin  control  system  for  Program  777.  In  this  configuration,  the 
spacecraft  will  be  operated  in  the  search  mode  with  additional  commands 
(i.e.,  counts)  sent  up  from  the  ground.  The  problem  with  operating  the 
spacecraft  in  the  normal  mode  is  that  when  the  platform  is  knocked  off 
the  earth  by  large  torque  disturbances,  the  on-board  logic  is  programmed 
so  that  after  3 invalid  signals,  the  controller  goes  into  standby.  To 
recover  the  spacecraft,  operator  intervention  involving  much  time  and 
effort  is  required. 

The  advantage  of  the  "fly-by-wire"  design  is  that  even  if  the  platform 
loses  earth  lock,  the  controller  will  continue  to  send  conmands  to  the 
motor  and  thus  reasonably  continuous  operation  of  the  spacecraft  can  be 
maintained.  The  use  of  the  narrow  coverage  antennas  is  precluded, 
although  the  1 to  2.5  degree  pointing  accuracy  achievable  by  the  "fly-by- 
wire"  technique  is  more  than  adequate  to  allow  the  use  of  the  earth 
coverage  antenna. 

Two  designs  for  the  ground  loop  controller  are  considered  - one  using 
proportional  control  only  and  one  using  proportional  plus  rate  control. 
Linear  analyses,  backed  up  by  nonlinear  simulation  results,  show  that 
the  best  performance  is  achieved  at  the  lower  RPM,  because  higher  position 
gains  can  be  tolerated.  In  effect  the  position  gains  selected  are  roughly 
1 count  per  1 degree  error  at  40  RPM,  and  1 count  per  2 1/2  degree  error 
at  75  RPM.  The  addition  of  the  rate  loop  produced  no  improvement  in  the 
steady-state  performance,  but  did  result  in  somewhat  smaller  transients 
in  response  to  large  disturbances  at  the  higher  RPM. 

It  is  thus  recommended  that  the  spacecraft  be  operated  at  the  lower 
RPM's.  In  addition  it  is  also  recommended  that  the  control  system  first 
be  tested  using  only  the  position  loop.  If  this  is  satisfactory,  the 
rate  term  need  not  be  included.  It  is  further  recommended  that  the 
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1 .0  INTRODUCTION 

The  777  spacecraft  #9433  despin  control  system  was  originally  designed 
to  point  the  platform  to  the  earth  with  an  accuracy  of  .15  degree.  In 
August,  1975,  large  pointing  errors  of  1.5  degrees  were  observed,  and  some- 
time later,  the  spacecraft  lost  lock  with  the  earth.  Since  then,  repeated 
attempts  have  been  made  to  bring  the  spacecraft  back  into  normal  mode,  but 
each  time,  sometimes  even  after  successful  operation  of  several  hours  or 
more,  the  spacecraft  is  knocked  off  the  earth,  and  normal  mode  control  is 
lost.  (The  problem  with  the  normal  mode  is  that  the  onboard  logic  is 
programmed  so  that  if  three  successive  platform  signals  are  received  out- 
side the  scanned  earth  chord,  the  controller  goes  into  standby.)  To  bring 
the  spacecraft  back  into  normal  mode  requires  operator  intervention,  and 
in  recent  months  this  activity  h»s  taken  a considerable  amount  of  time. 

The  cause  of  this  anomalous  performance  has  been  the  subject  of  much 
study  in  the  past  few  months  (Reference  1).  Basically  the  anomaly  is 
characterized  by  an  increase  in  the  constant  running  friction  {-  100  in-oz) 
across  the  rotor-platform  interfaces,  accompanied  by  occasional  torque 
transients.  It  was  concluded  that  the  anomaly  was  most  likely  due  to  poor 
lubrication  in  the  DMA,  although  the  lubrication  supply  system  is  not  very 
well  understood. 

It  has  recently  been  suggested  that  an  attempt  be  made  to  operate  the 
spacecraft  in  the  rate/search  mode,  with  a ground  commandable  loop  to  the 
spacecraft.  In  doing  so,  it  was  realized  that  the  pointing  accuracy  would 
be  degraded,  but  the  advantage  in  this  type  of  operation  is  that  even  if 
the  platform  is  knocked  off  the  earth,  the  controller  can  be  progranmed  so 
that  it  will  not  automatically  "quit”  as  in  the  case  of  the  normal  mode, 
but  will  continue  issuing  commands  to  the  i.'Kitor.  Thus,  reasonably  continuous 
operation  of  the  spacecraft  can  be  maintained. 

This  document  examines  the  possibilities  of  the  ground  loop  in  the 
control  system,  and  suggests  a control  system  that  can  keep  the  steady-state 
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pointing  accuracy  to  within  1 to  2 1/2  degrees.  There  are  three  main 
obstacles  to  finer  pointing  accuracy  by  using  the  rate/search  mode  and  a 
ground  loop.  The  first  is  the  rate  quantization  already  existing  in  the 
rate  loop  in  the  control  system  of  the  spacecraft.  The  other  two  are  con- 
cerned with  the  ground  loop  - first,  the  transport  lag  which  is  currently 
estimated  at  between  2-6  seconds,  and  second,  there  is  quantization  in  the 
ground  loop  as  well.  That  is,  only  an  integer  number  of  counts  can  be 
sent  from  the  ground  to  the  spacecraft,  and  this  quantization  turns  out 
to  be  quite  large. 

In  order  to  fully  understand  the  effects  of  quantization,  a brief 
look  is  taken  at  the  control  system  with  no  delay  in  the  ground  loop. 

This  is  done  in  Section  2.0,  and  the  ground  loop  control  system  consists 
simply  of  proportional  control.  The  effect  of  varying  the  spin  speed 
(40  to  75  RPM)  is  also  considered.  Then  in  Section  3.0,  the  effect  of 
the  transport  lag  in  the  ground  position  loop  is  studied.  Finally,  in 
Section  4.0,  a ground  loop  consisting  of  proportional  plus  rate  control 
is  examined.  Implementation  considerations  of  the  ground  loop  are  discussed 
in  Section  5.0.  Conclusions  and  Recommendations  are  discussed  in  Sec- 
tion 6.0.  Appendix  A is  a summary  of  the  results  of  an  analysis  of  the 
time  delays  to  be  expected  in  the  ground  loop,  and  Appendix  B describes 
the  simulation  used  in  this  document. 
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2.0  A "ZERQ-PELAY“  SYSTEM 

This  section  examines  the  control  system  in  the  rate/search  mode  with 
a ground  commandable  loop  which  consists  only  of  a gain,  Kp,  times  position. 
The  2-plane  transfer  function  will  be  determined  in  order  to  compute  the 
range  of  Kp  for  stable  operation.  Then,  some  non-linearities  of  the  system 
are  examined,  with  time  traces  plotted  for  various  RPM.  All  this  will  be 
done  assuming  zero-delay  in  both  the  rate  loop  and  the  position  loop. 


2.1  Z-Plane  Analysis  of  the  Linear  System 

A linear  model  of  the  despin  controller  in  the  rate/search  mode  with 
the  ground  loop  in  is  shown  in  Figure  2-1.  The  open  loop  transfer  func'. 'on 
L{s)  assuming  zero  transport  lag  is 


L(s) 


^(l-e  * ) 


kc*(s  + 


B) 


, ''s  4 ''c'  "m  {s  + B)  (s  + Kp/(Kj^T2)) 

T^  d-e’^^)  (s  + B)  (s  + A) 


where 


A ■ Kp/(K|,T2) 

T * spin  period  (»  1 sec  at  60  RPM) 
Tj  » 2T  ■ sampling  period 


(2-1) 
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Table  2-1.  Parameter  Values  for  Figure  2-1 


Variable 


Value 

Description 

.0667 

Filter  breakpoint  (rad/sec) 

72 

Platform  inertia  (slug-ft  ) 

.00586 

D/A  converter  gain  (volts/count) 

G.S.  1 1.00 

G.S.  2 1.39 

G.S.  3 2.01 

G.S.  4 2.82 

G.S.  5 3.98 

G.S.  6 5.65 

G.S.  7 7.97 

G.S.  8 11.32 

Selectable  gain  (volts/volt) 

.4 

Motor  gain  (ft-lb/volt) 

5.51 

Compensation  gain  (volt/volt) 

40  RPM  - 1.5 

60  RPM  - 1.0 

75  RPM  - .8 

Spin  period  (sec) 

2T 

Sampling  period 

See  Appendix  A 
(2-6  sec) 

Delay  in  ground  loop 
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The  Z-transform  of  L(s)  is  determined  by  partial  fraction  expansion. 
Without  the  measurement  delay,  the  controllers'  open  loop  Z-plane  transfer 
function  is 


Uz)  = Z [L(s)] 


K — Z 
0 z 


~(s  B)  (s  + aY 


K„  ii 
0 z 


= K ^ 
■^0  ' z 


rxi  + C2  + C3  + C4-1 

s 7 7 7 

(2-2) 

c,  z + C2  TgZ  + C3  J^^ziz  + 1)  + C4  T^^z  (z^  + 4 z + 1)-] 

T 


z-1  (z-1)^  2(z-l)^ 


6(z-l)‘ 


where 


so 


C4  = AB 
03=  A + B 

C2  = 1 

Cl  = 0 


L(z)  = <0  ^ 


6T  z(2^  - 2z+l)  + (3A  + 3B)T  ^z(z^-l)  + ABT  ^z(z^  + Az+l ) 


6(z-l)‘ 


^0 


z''(6Tg  + (3A  + 3B)  + ABT^)  + z(-12Tg  + 4 ABT  + 6T  -(3A+3B)T  ^+ABT 


-^+ABT  ^ 


6(z-l)' 


L J 


(2-3) 


where  = root  locus  gain  in  the  Z-plane. 
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The  two  zeros  of  the  Z-plane  transfer  function  are  B and  a.  The  zero 

BT 

B which  is  associated  with  the  S-plane  filter  breakpoint  B,  (b  = e"  s)  is 
independent  of  A.  The  second  Z-plane  zero,  a,  is  plotted  in  Figure  2-2  as 
a function  of  A (at  60  RPM,  or  T=l).  The  zero  a approaches  e"^^^  as  A becomes 
smal  1 . 


The  Z-plane  root  locus  associated  with  the  transfer  function  at  60  RPM 
(T=l)  is  shown  in  Figure  2-3.  Very  roughly,  the  root  locus  (the  portion  off 
the  axis)  is  a circle  with  a radius  equal  to  1-a.  It  is  thus  evident  that  in 

order  for  most  of  the  root  locus  to  be  inside  the  unit  circle,  that  o be  > 0. 
From  Figure  2-2  this  implies  that  A be  <1  (for  60  RPM).* 

Similar  plots  of  a as  a function  of  A were  made  at  75  RPM  (T=.8,  or  T^  = 
1.6)  and  40  RPM  (T  = 1.5,  T^  ="3).  These  are  shown  in  Figure  2-4.  These 
show  that  at  the  lower  RPM's  the  s-plane  zero  A must  be  reduced. 


It  is  worthwhile  at  this  point  to  digress  for  just  a moment  to  consider 
what  the  value  of  A (the  S-plane  zero)  really  means,  in  terms  of  the 
ground  commandable  loop.  The  gain  1(  in  this  position  loop  is 


K 


P 


(2-4) 


2 

which  means  that  Kp  counts  should  be  sent  for  every  AKj^T 
Since  = 261,  this  means  that 


radians  of  error. 


Kp  = 261  A T^  (2-5) 

2 

or  in  terms  of  degrees,  4.56  AT  counts  should  be  sent  for  every  degree 
of  error.  At  60  RPM  {T=l)  typical  values  of  A turn  out  to  be  on  the  order  of 
.1  or  .2.  A value  of  A = .1  means  .456  counts  for  every  degree  of  error,  or 
since  only  an  integer  number  of  counts  can  be  sent,  _]_  count  should  be  sent 


* This  turns  out  to  be  a large  upper  bound;  in  fact  values  of  A = .1  and  .2 
were  used. 


EMI-LOGARITHMIC  46  5BI3 


SIMI  LOGARITHMIC  46  5813 
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for  every  2 degrees  of  error.  It  has  been  shown  that  at  the  lower  RPM,  the 

value  of  A decreases,  but  since  T increases,  the  effect  on  K is  not  obvious 

P 

at  this  point.  Clearly,  however,  the  accuracy  to  be  expected  from  this  type 
of  system  is  going  to  be  on  the  order  of  degrees,  and  not  fractions  of  degrees. 

A gain  phase  plot  of  the  system  with  A = .1  and  A = .5  is  shown  in 
Figures  2-5  and  2-6  for  60  RPM  (0  delay).  These  two  figures  show  that,  at  60  RPM 
a higher  value  of  A is  possible  if  the  gain  state  is  increased.  In  other  words, 

A is  dependent  on  the  gain  state  chosen,  and  in  order  to  calculate  the  "best 
A",  one  should  choose  the  optimum  gain  state  for  each  RPM.  This  suggests 
formulating  a new  problem,  where  A is  the  overall  gain  in  the  characteristic 
equation.  Thus  a root  locus  will  be  plotted  where  the  parameter  A is  the 
variable  gain.  The  approach  used  will  be  to  rewrite  the  characteristic  equation 

1 + L^(2)  + A Lp(z)  = 0 

where 

Lp(z)  is  the  rate  loop  transfer  function 
and  ^Lp(z)  is  the  position  loop  transfer  function 

A L (z) 

as  1 + ^ = 0 

1 » L„(Z) 

SO  that  the  new  "open  loop"  transfer  function  to  be  plotted  is  H(z) 


(2-6) 


(2-7) 


With  no  delay,  the  transfer  function  of  the  inner  rate  loop  is 

.2 


L„(z)  ■ V 


o 


z-1 

•T 

1 . B 

z 

• L 

^ IT 

2-1 

Z 


T 2 B tS  (z  + 1) 

-5—-  + ^ — 

/ _ 1 \0 


L (i-l)'^  2 (z  - 1)' 


2 2Tj(2-1)  + BT^(z+l) 

2(z-1)2 


^ (BT^  > 2T3)  - 2T^  . BT^ 


wnere 


02  = BT;.2T^ 


Oi  = bt;  - 21^ 
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Lr(2)  . KT 


2 0 2 ^ <3^ 

~2(z-1)2 


(2-8) 


The  Z- transform  of  the  position  loop  is 
Al-n(2)  = A ^ Z 


0 z 


= KT^A  lill  (z.l)  . B T^zU^  * 4z*i; 

° ^ L 2(z-i)^  eYz-i)'^ 


T^A  ^^s(z^-l)  •*•  BTs(z^+4z+l 


6(z-l)' 


ALp(z)  = KJ^A  ^2  ^ ^ *‘1  z ^ 

L 6 (z-1)^ 


(2-9) 
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where 

3T^  + BT^ 
s s 

4BT^ 

BT^  - 31^ 

The  characteristic  equation  is 

1 + Lj^(z)  ^ ALp{z)  = 0 (2-10) 


or  substituting  in  the  expressions  for  Ln(z)  and  L (z), 

K p 


1 + K, 

t2 

a^Z  + o-j 

+ K T^A 

( 

D 

[ 2(z-1)2_ 

0 

Multiplying  through  by  6(z-l)^ 

gives 

6(z-l) 

3(z-l) 

• [a. 

2Z  + a^]  + 

KqT^  a [p 

This  must 

be  rewritten  in 

the  form  of  l+AH(z)  = 0 

so 

1 + 

A 

hr^z 

2 ^ 

+ W-jZ  + 

“ol 

(Z-1)  [6(z-1)2  + 

3 j . 

*■  O])] 

1 + 

A i: 

0 

— — — * 

■ + z + 

^2  z u^z  + u 


6 (z-1)' 


= 0 (2-11) 


(z-1)  [6  + (3K  a„  - 12)  z + 3 K a,  + 


2 " ‘ ^ J ♦ 6] 


= 0 


(2-13) 


It  is  clear  from  this  trahsfer  function,  that  the  value  of  A depends  on  K 
(i.e.,  the  gain  state)  and  T^. 
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It  will  be  useful  at  this  point  to  reintroduce  the  parameter  Kp,  which 
is  the  ground  position  loop  gain,  since  the  parameter  A is  somewhat  arti- 
ficial. Recall  that  A was  the  S-plane  zero  in  the  zero  delay  system; 
however,  when  the  characteristic  equation  is  written  as  above,  it  makes 
little  sense  to  isolate  A,  since  the  position  gain  Kp  is  dependent  on  the 
product  AT^.*  Let  the  variable  Kp'  be  defined  as 

= VSi  • 

2 

and  substitute  K ' for  AT  in  the  above  equation.  Then  the  characteristic 
P 

equation  becomes 


(z-1)  [6  z^  + (3KJ^  02-12)  z + 3 a.,+6] 


(2-15) 


Figure  2-7  contains  the  root  locus  of  H(z)  for  60  RPM  (T»l)  for  gain 

states  5-8.  The  zeros  are  independent  of  the  gain  state  chosen.  One  of 

the  zeros  is  the  zero  associated  with  the  S-plane  break  frequency  B (i.e., 
-BT 

z = e s)  while  the  other  is  located  in  the  left  half  plane  just  outside 

the  unit  circle  on  the  real  axis.  There  are  3 Z-plane  poles  - one  at  z=l , 

and  two  dependent  on  the  gain  state  chosen.  One  of  these  latter  poles 

BT 

approaches  the  Z-plane  zero  at  e while  the  other  moves  to  the  left 
with  increasing  gain  state.  As  Figure  2-7  shows,  a higher  value  of  Kp' 
is  possible  at  the  higher  gain  states.  However,  the  association  of  a 
higher  value  of  Kp'  with  a higher  gain  state  is  true  only  up  to  a point. 
Figure  2-8  shows  the  root  loci  at  40  RPM  for  various  gain  states,  and  as 
can  be  seen  in  gain  state  7 (and  also  in  8),  the  system  is  unstable  for  all 
values  of  Kp'.  This  is  because  the  Z-plane  pole  has  moved  so  far  to  the 
left  that  it  has  passed  the  Z-plane  zero  outside  the  unit  circle. 


0 

f 

r. 

w 

' « 
f 


* In  addition,  in  the  next  section,  when  the  system  is  considered  with 
different  delays  in  each  loop,  the  parameter  A loses  its  significance 
as  the  s-plane  zero,  and  it  retention  only  leads  to  confusion. 
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Transfer  Function  H(z) 


poles  - (.80.  + .10),  (1.0) 
GAIN  STATE  5 


poles  - (.81,  0),  (.62,  0),  (1,  0) 
GAIN  STATE  6 


poles  - (.84,  0),  (.36,  0),  (1,0)  poles  - (.86,  0),  (.01,  0),  (1,  0) 

GAIN  STATE  7 GAIN  STATE  8 


* Zeros  at  (.88,  0),  (-1.05,  0) 


Figure  2-7.  Root  Locus  of  H(z)  as  a Function  of 
Gain  State.  60  RPM 
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Figure  2-9  shows  the  maximum  values  of  K„'  = K_/K„  plotted  as  a function 

P P ^ 

of  the  gain  states,  for  various  RPM.  This  answers  one  of  the  questions 
posed  earlier  in  the  section  regarding  the  position  gain  Kp  = 261  Kp'.  The 
gain  Kp  meant  that  4.56  Kp’  counts  should  be  sent  for  every  degree  of  error, 
or  conversely,  since  only  an  integer  number  of  counts  can  be  sent,  1 count 
should  be  sent  for  every  1/4.56  Kp'  degrees  of  pointing  error.  Figure  2-9 
shows  that  at  the  lower  RPM,  a higher  value  of  Kp'  is  possible.  This 
suggests  that  finer  control  can  be  obtained  at  the  lower  RPM's.  However, 
it  should  also  be  noted  that  at  the  lower  RPM's,  the  control  signal  is 
sent  less  often.  These  effects  will  be  tested  by  simulation  at  the  end  of 
this  section. 

Figure  2-9  also  shows  that  at  60  and  75  RPM  the  spacecraft  should  be 
operated  in  the  rate/search  mode  at  gain  states  7/8  while  using  the  ground 
loop,  and  at  40  RPM,  somewhat  lower  gain  states  (perhaps  5 or  6)  should  be 
used.  More  will  be  said  about  this  when  the  gain  phase  plots  are  drawn. 

This  linear  analysis  could  be  carried  out  even  further  to  determine  the 
ideal  damping  ratio  and  its  dependence  on  A,  etc.,  but  before  carrying  on 
an  extensive  analysis,  it  is  worthwhile  to  examine  some  of  the  non-linearities 
in  the  system  to  see  how  they  affect  the  system's  performance. 

2.2  Non-Linearities  in  777  Simulation 

Figure  2-10  is  a diagram  of  the  nonlinear  control  system.  The  signifi- 
cant nonlinearities  are  discussed  below  to  determine  their  effect  on  the 
systems  linear  model. 

2.2.1  Counter  Quantization  - Rate  Loop 

In  the  rate  loop,  one  count  is  sent  for  every  .22/T^  deg/sec  rate  error. 
This  means  that  the  effective  rate  quantization  is  as  follows: 
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RPM 

T 

75  RPM 

0.8 

1 count  = .22/(.80)^  = .34  deg/sec 

60  RPM 

1.0 

1 count  = .22/1  = .22  deg/sec 

40  RPM 

1.5 

1 count  = .22/(1. 5)^  = .098  deg/sec 

With  this  large  rate  quantization,  it  means  that  it  will  be  impossible 
for  the  pointing  error  of  the  spacecraft  to  settle  down  and  converge  as  in 
the  normal  mode.  What  can  be  hoped  for,  however,  is  that  the  steady-state 
error  can  be  kept  within  a limit  cycle  on  the  order  of  1 or  2 degrees. 

The  above  table  also  shows  that  at  the  lower  RPM,  since  the  rate  quantiza- 
tion is  smaller,  finer  pointing  error  control  can  be  maintained. 

2.2.2  Position  Loop  Quantization 

The  other  significant  nonlinearity  in  the  control  system  is  the  quanti- 
zation in  the  position  loop  which  has  already  been  discussed  in  Section  2-1. 
Only  an  integer  number  N of  counts  can  be  sent  from  the  ground  computer, 
and  this  number  is 

N = 4.56  Kp'  counts/degree  of  error 

or  similarly 

1 count  for  every  1/4.56  Kp'  degrees  of  error. 

It  is  interesting  to  study  the  effects  of  these  nonlinearities  on  the 
control  system.  Figure  2-11  is  a plot  at  60  RPM  of  the  response  to  a torque 
disturbance  of  10  in-oz  for  20  sec.  (There  was  no  delay  in  either  loop  in 
this  test  case.)  With  no  position  loop  quantization  or  rate  quantization, 
the  error  settles  down  to  zero  as  in  the  normal  mode.  With  the  position 
quantization  in  (but  still  no  rate  quantization)  there  is  a steady-state 
error  no  larger  than  1/4.56  Kp'  degrees.  However,  with  no  position 
quantization,  but  with  rate  quantization,  we  see  the  oscillatory  response 
previously  predicted.  What  this  means  is,  with  the  present  spacecraft  even 
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if  the  ground  position  loop  (consisting  of  a gain)  were  ideal,  i.e.,  no 
delay  and  no  quantization,  it  would  still  be  impossible  to  get  the  steady- 
state  error  to  settle  down. 

Figure  2-12  shows  the  pointing  error  response  with  all  quantization 
in  (i.e.,  both  rate  loop  and  ground  position  loop)  but  still  no  delay  in 
either  loop.  After  an  initial  transient,  the  error  oscillates  between 
+ 1/4.56  Kp' , or  + 1.1  degrees  for  the  case  of  Kp'  = .2.  Also  shown  in 
this  figure  are  the  ground  loop  commands,  the  rate, and  the  rate/error  com- 
mands due  to  the  rate  loop  in  the  spacecraft. 

In  Figures  2-13  and  2-14,  the  effect  of  increasing  the  position  gain 
is  examined,  and  as  predicted,  the  limit  cycle  amplitude  decreases  with 
increasing  gain,  so  that  in  Figure  2-14,  the  limit  cycle  is  .43  degree  at 
60  RPM,  for  Kp'  = .5  (zero-delay  system). 

Consider  now  the  case  of  40  RPM,  where  the  rate  quantization  in  the 
spacecraft  is  smaller  (1  count  per  .10  deg/sec  rate  error).  Figure  2-15 
shows  the  effects  of  removing  quantization  from  both  loops  and  then  each 
loop  separately,  and  the  effect  of  the  smaller  rate  quantization  can  be 
seen  in  the  lower  plot  in  Figure  2-15,  where  the  limit  cycle  is  somewhat 
smaller  than  the  corresponding  case  at  60  RPM  (Figure  2-11).  In 
Figure  2-16,  with  quantization  in  both  loops,  the  pointing  error 
oscillates  between  + 1/4.56  Kp'  (+  1 deg  for  Kp'  = .2)  as  in  the  case  at 
60  RPM.  The  effect  of  increasing  the  gain,  in  the  lower  limit  cycle,  is 
again  shown  in  Figure  2-17. 

At  75  RPM,  the  effects  of  the  larger  rate  quantization  are  clearly 
seen  in  the  third  plot  of  Figure  2-18,  where  the  position  quantization 
has  been  removed.  Typical  responses  of  the  system  showing  the  effects 
of  increasing  the  gain  on  the  limit  cycle  behavior  are  illustrated  in 
Figures  2-19  and  2-20. 

The  purpose  of  this  section  on  the  zero-delay  system  has  been  to 
illustrate  the  effects  of  the  rate  and  position  quantization  on  the 
system,  and  show  that  even  if  the  ground  loop  were  ideal,  i.e.,  no  delay, 
it  would  be  impossible  to  get  a convergent  response  of  the  pointing 


BD.nn  (ED.nn  HHO.nn  ian.nn  hdd.qd 

Time  (sec) 


76.7341.9-027 
Page  36 


error  by  using  the  search  mode  and  a ground  controller.  The  steady-state 
error  response  is  always  going  to  be  oscillatory,  and  the  next  section 
will  show  that,  when  the  effects  of  the  delay  are  considered,  the  minimum 
limit  cycle  to  be  expected  will  be  about  1 deg.  It  will  also  show  that 
the  lower  RPM's  (a  range  of  40  to  75  was  considered)  provide  somewhat 
better  limit  cycle  response  due  to  the  fact  that  a higher  value  of  the 
ground  position  gain  is  possible  at  the  lower  RPM's,  and  also  somewhat 
due  to  the  decreased  rate  quantization  at  the  lower  RPM's.  (The  former 
fact  was  indicated  by  Figure  2-9.) 


f 

f 
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3.0  DESIGN  OF  A POSITION  LOOP,  CONSIDERING  THE  EFFECTS  OF  THE  DELAY 

In  this  section,  the  analysis  of  Section  2 will  be  repeated,  but  the 
delay  in  both  the  rate  loop  (=  .5T)  and  the  position  loop  will  be  considered. 
Currently,  the  position  loop  delay  is  estimated  to  be  about  2-6  sec., 
depending  on  the  spin  speed.  (See  Appendix  A.) 

Before  going  into  the  analyses,  the  main  conclusions  from  Section  2 
(zero-delay  system)  will  be  repeated  because  they  provide  guidelines  for 
the  type  of  behavior  to  be  expected  from  the  system  with  transport  lag 
included.  The  main  conclusions  v<ere; 

• Rate  quantization  is  the  main  obstacle  to  obtaining  a convergent 
response  of  the  pointing  error,  and  instead  leads  to  an  oscillatory 
response. 

• There  is  less  rate  quantization  at  lower  RPM's  so  possibly  better 
performance  at  the  lower  spin  speeds  can  be  expected. 

• The  amplitude  of  the  limit  cycles  of  the  zero  delay  system,  taking 
into  account  quantization,  can  be  predicted  by  1/(4. 56  Kp* ) degrees 
(assuming  stable  value  for  Kp'.) 

• The  linear  stability  analysis  results  show  that  higher  values  of 
the  position  gain  Kp  = 4.56  Kp'  are  possible  at  lower  RPM's,  so 
so  perhaps  a smaller. limit  cycle  can  be  obtained. 
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3.1  Linear  Analysis 


The  open  loop  Z-plane  transfer  function  of  the  controller's  inner  rate 
loop  is 

-sT, 


= K t2  .(z-1J  2 e 

0 z 


lllil 

,3 


(3-1) 


The  Z-plane  transfer  function  is  evaluated  by  partial  fraction  expansion  and 
modified  Z transforms. 


Lr  (z) 


= K T' 


l-l 

^ J 

¥)&  • 

^s 

(z-l)2 

(^2  * 

j and  m = .75 

J 

O T 2 


'1  's  m' 


m . 2ni+1  . 2 

^ (z-l)2  (z-1)^ 


2z(z-l)^ 


(3-2) 


where 


Tg  (2m  + BT^m'") 


Tj  [2-4m  - 2BT5m'^  + BT^(2m  + 1)] 


Tj  [2m-2  + BTj(m'^  - 2iiH-l)] 


= [2(m-l)  + BT5(m-l)''] 


76.7341.9-027 
Page  39 


The  Z-transform  of  the  position  loop  is 


Lp(2)  • KpKp’ 


(3-3) 


Because  the  delay  in  this  position  loop  is  on  the  order  of  several 
seconds,  i.e.,  large  compared  to  the  sampling  interval,  T^,  it  is  convenient 
to  study  cases  where  the  delay  is  a multiple  of  the  sampling  period.  If 
this  is  done,  then  the  Z-transform,  L (z)  is  easily  taken  by  multiplying 
the  z transform  of  Lp(z)  with  no  delay  by  1/z  , i.e.. 


Lp<^)  ■ Vp'  ^ 


if) 


(3-4) 


where  k is  the  number  of  sampling  periods  of  the  delay  in  the  ground  loop. 

The  transform  of  the  position  loop  without  delay  has  already  been  evaluated  in 
Section  2.  Therefore,  the  transfer  function  of  the  position  loop  with  delay 


^ >^1  ^ ^ % 


6 (z-1) 


(3-5) 


where 


3T^  + BT^ 
s s 

4BT^ 

s 

BT^  - 3T^ 
s s 


The  transfer  function  L(z)  of  the  combined  system  is 
L(z)  = Kp'  Lp(z)  + Lr(z), 


K ' K /u«  Z^  + u,  Z + y \ 

L(z)  - ^ 

\ 6(z  - 1)^  / 


2 '"2 


(no  z'^  + ni  z + n„) 


2z(z  - D' 


(3-6) 
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The  poles  of  the  combined  system  K ' L (z)  + Ld(z)  are: 

P P “ 


1.  three  poles  at  z = 1 

2.  1 pole  at  origin  for  k = 0,  1 

3.  k poles  at  origin  for  k = 2, 


The  zeros  of  the  combined  system  are  not  as  easily  solved  for.  One  of 
the  zeros  is  associated  with  the  S-plane  filter  breakpoint  B,  or  z = e 
The  location  and  number  of  the  other  Z-plane  zeros  is  dependent  on  the  ground 
loop  delay,  and  the  ground  loop  position  gain  K^'. 

In  order  to  solve  for  the  value  of  K^'  to  use,  which  depends  on  the 
gain  state  chosen  and  the  spin  speed,  the  characteristic  equation  of 
the  system  will  be  rewritten  so  as  to  isolate  Kp'  as  the  overall  system 
gain,  as  was  done  in  Section  2.0.  That  is,  if  Lp(z)  is  the  transfer 
function  of  the  inner  rate  loop,  and  Kp'(z)  is  the  transfer  function  of 
the  position  loop,  the  characteristic  equation 


will  be  rewritten  as 


1 + 4(z)  + Kp'(z)  = 0 


, I'd'  I-d*") 


(3-7) 


(3-8) 


and  the  new  open  loop  transfer  function,  H(z) 


K'  >-„(2) 

■ A 

will  be  plotted  in  the  Z-plane. 

The  transfer  function  of  the  inner  loop  is 

2z(z-1)2 


(3-9) 


(3-10) 
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and  the  transfer  function  of  the  outer  loop  is 


K • L (z)  = 

P P 2^ 


*^2  ^ ^ h ^ ^ ^0 
6(2-1)^ 


(3-11) 


When  the  characteristic  equation  1 + Ld(z)+  K ‘ L (z)  = 0 is 

K p P 

rewritten  as 


KJ  k(2) 

TI7 


(3-12) 


the  result  is 


1 + 


z'^“^z-l)  [6z^  + z^3y^2  - ■'2)  + z(6  + 3KJ^n^)  + J J 


= 0 


or 


1 + G(z)  = 0 (3-13) 

where  G(z)  is  equal  to  the  expression  in  brackets  in  the  previous  equation. 

The  root  locus  of  G(z)  was  considered  for  different  gain  states  at 
each  RPM  with  various  delays  in  the  ground  loop  (i.e.,  different  values 
of  k),  where  the  ground  loop  delay  was  considered  to  be  multiples  of  the 
sampling  period.  The  maximum  values  of  the  previous  gain  Kp'  for  stability, 
i.e.,  those  values  which  cause  the  roots  of  the  characteristic  equation  to 
be  on  the  unit  circle,  are  shown  in  Figure  3-1. 
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3.2  Position  Loop  Results 

The  range  of  the  position  loop  gain  Kp  for  stability  of  the  non-linear 
system  is  very  nearly  predicted  by  the  stability  of  the  linear  system  as  shown 
in  Figure  3-1  if  the  position  loop  counter  is  implemented  as  follows  (Figure  3-2): 


Counts 


1 

-1 

1 K 0 

P 

-1 

Figure  3-2.  Implementation  of  Ground  Loop  Counter 

The  region  around  the  origin,  i.e.,  "0  counts"  is  of  particular  importance  because 
it  has  already  been  shown  in  Section  2.0  that  the  steady  state  error  exhibits  a limit 
cycle  behavior  oscillating  between  + 1/Kp  degrees,  with  the  position  loop  COTimanding 
at  most  one  count  in  either  direction.  (The  results  for  Section  2.0  considered 
a zero-delay  system,  but  the  limit  cycle  of  the  system  with  the  ground  loop  delay 
is  only  slightly  larger  than  1/Kp  degrees,  and  still  only  one  count  is  needed  in 
the  steady  state.)  It  turns  out  that  it  is  possible  for  the  control  system  to 
cotimand  more  than  one  count  in  response  to  certain  transients;  however,  even  then, 
there  are  certain  advantages  to  limiting  the  counter  to  ^ 1 count,  and  these  will 
be  presented  later  in  this  section.  The  important  point  is  that  any  gain  selected 
should  be  done  so  with  the  ground  commanded  counter  implemented  as  in  Figure  3-2. 

Simulation  results  indicate  that  because  the  effective  gain  is  somewhat 
reduced  by  the  quantization,  somewhat  higher  values  for  the  gains  could  be 
tolerated  than  those  indicated  as  being  the  maximum  value  in  Figure  3-1,  especially 
at  the  lower  RPM.  The  following  values  of  Kp'  were  chosen  with  the  idea  of 
i picking  a higher  value  of  K ' in  order  to  reduce  the  limit  cycle. 
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40 

RPM 

■ 

.3 

60 

RPM 

K ' = 

P 

.15 

75 

RPM 

K ' = 

.10 

P 

The  gain  phase  plots  of  the  transfer  function  of  the  control  system,  whose  open 
loop  transfer  function  is  L(z), 


where 


L{z)  = ^ 


z(e 


-T5S/4 


(s+B) 


and  ^s’  ^ ~ 

are  shown  in  Figures  3-3  through  3-10  for  various  ground  loop  delays 

at  each  RPM.  The  final  selection  of  gains  in  terms  of  counts  per  degree,  and 

also  the  ideal  gain  state  at  each  RPM  is  shown  in  Table  3-1. 


Table  3-1.  Position  Loop  Gain  Selection 


RPM 

GAIN  STATE 

GAIN  (COUNTS  PER  DEGREE) 

•40  RPM 

6 

1.368 

60  RPM 

8 

.684 

75  RPM 

8 

.456 

The  counter  should  be  limited  to  +1  count,  with  1 count  being 

sent  when  the  magnitude  of  KpS  exceeds 

1. 

DELAY  = 1.6  SEC 


Figure  3-8.  777  Search  Mode  with  Ground  Position  Loop 


DELAY  = 3.2  SEC 


Figure  3-9.  777  Search  Mode  with  Ground  Position  Loop  75  RPM 


3.3  Performance  Analysis 


In  analyzing  the  performance  of  the  control  system,  two  criteria  were 
examined  --  the  peak  pointing  error  during  a transient,  and  the  limit  cycle 
performance.  At  each  RPM,  the  worst  case  delays  from  Appendix  A were  assumed. 

3.3.1  Steady-State  Performance 

Figure  3-11  shows  the  response  of  the  system  at  40  RPM  using  the  gains  and 
gain  state  from  Table  3-1  to  a simulated  disturbance  of  10  in-oz  for  20  sec. 

The  limit  cycle  is  + .8°,  with  the  ground  loop  coimanding  at  most  + 1 count. 
(There  was  actually  a limit  on  the  counter  of  5 counts  in  this  case.) 

Figure  3-12  shows  the  response  of  the  control  system  to  the  same  disturbance 
at  60  RPM,  and  a limit  cycle  of  under  2°  is  obtained. 

Although  in  both  of  the  previous  cases  a true  oscillatory  response  was 
obtained  in  the  steady-state,  i.e.,  the  control  system  went  from  one  side  of 
the  deadzone  to  the  other,  it  is  also  possible  for  the  system  to  stay  on  only 
one  side  of  the  deadband.  This  is  shown  in  Figure  3-13,  where  a different 
initial  condition  was  used  at  40  RPM.  In  this  case,  the  ground  loop  is  com- 
manding counts  on  only  one  side,  after  the  initial  transient.  (The  rate  loop 
in  the  spacecraft  is  also  issuing  counts.)  Still,  the  maximum  peaks  on  the 
one  side  are  at  .8°,  as  before  (40  RPM),  and  this  behavior  is  what  is  referred 
to  as  "limit  cycle"  in  this  report. 

3.3.2  Effect  of  Limiting  Counter 

A large  disturbance  (50  in  oz  for  20  sec)  was  simulated  in  Figures  3-14a  and 
3-1 4b  (60  RPM).  In  Figure  3- 14a,  the  counter  was  limited  to  + 5 counts  (only  + 

4 counts  were  actually  commanded)  and  in  Figure  3-14b,  the  counter  was  limited 
to  + 1 count.  The  important  point  to  notice  is  that  the  initial  peak  pointing 
error  was  the  same  in  both  cases  (»4.8°)  but  that  in  the  case  where  the  counter 
had  limits  of  + 5,  a large  overshoot  was  seen  (*^6°)  whereas  in  the  case  where 
the  counter  was  limited,  the  overshoot  was  greatly  reduced.  The  large  overshoot 
in  the  first  case  was  due  to  the  time  delay  in  the  system.  It  seems  reasonable, 
then  to  limit  the  counter  to  + 1 count,  in  order  to  eliminate  these  large  over- 
shoots especially  in  light  of  the  fact  that  in  the  steady  state,  1 count  is  the 


1 count  per 
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most  needed  anyway.  Another  advantage  is  that  in  limiting  the  counter,  if 
the  gain  selected  is  too  high,  the  system  will  not  go  unstable,  but  merely 
limit  cycle  at  a higher  frequency  and  somewhat  higher  amplitude  than  if  a 
stable  gain  had  been  chosen.  An  example  of  this  is  provided  in  the  next 
three  figures,  which  show  a test  case  at  75  RPM. 

Figure  3-15  shows  the  response  of  the  system  to  the  same  disturbance  used 
earlier  (10  in-oz  for  20  sec)  and  with  a stable  value  of  Kp  selected, 

(Kp'  = .1  implies  1 count  for  every  2.2°)  and  steady-state  peaks  of  about 
2.4°  are  obtained,  and  the  counter  commands  at  most  1 count.  (Note  in  this 
case  that  the  responses  is  again  "one-sided".)  In  Figure  3-16  an  unstable 
value  of  Kp'  was  purposely  selected  (Kp'  = .2,  or  1 count  for  every  1.1°  error) 
and  the  counter  limits  were  placed  at  + 5.  The  system  response  is  obviously 
unstable,  and  limit  cycles  at  + 16°  due  to  the  counter  limit.  However,  if  the 
counter  limits  had  been  placed  at  + 1 , as  shown  in  Figure  3-17,  the  error 
oscillates  between  ^3°.  Obviously,  this  is  not  optimum  either,  because 

^ better  response  was  obtained  in  Figure  3-15,  where  a 1 ower , stable  variable 

of  gain  was  selected,  but  what  should  be  noted  is  that  limiting  the  counter 
sort  of  "desensitizes"  the  control  system  and  acts  as  an  insurance  policy 
against  the  consequences  of  too  high  a gain  selection. 

3.3.3  Effect  of  Rate  Loop  Offset 

One  other  effect  that  can  change  the  character  of  the  steady-state  per- 
formance is  if  the  counter  in  the  rate  loop  is  not  "dead  set"  around  zero. 

The  ideal  model  for  the  rate  loop  counter  is  shown  below  in  Figure  3-18a, 
where  the  point  of  zero  rate  is  in  the  center  of  the  0 count.  However,  if 
the  rate  loop  counter  is  modeled  as  in  Figure  3-18b,  this  changes  both  the 
amplitude  and  frequency  of  the  pointing  error  response.  This  is  shown  in 
Figure  3-19,  where  an  offset  of  .3  was  used  instead  of  .5  at  75  RPM.  Also 
shown  in  Figure  3-19  are  the  cormiands  issued  by  the  rate  loop  in  the  space- 
craft. 

i 


Figure  3-17  Response  of  Control  System  at  75  RPM  [:  ;h 

Unstable  Gain  Selected,  Counter  Limits  at  + 1 


0 


76.7341.9-027 
Page  64 


Figure  3-18a.  Ideal  Rate  Counter  Figure  3-18b.  Rate  Counter  with 

(Offset  = .5)  Offset 

3.3.4  Transient  Response 

The  second  performance  criteria  was  the  peak  transient  response  to  sudden 
disturbances.  Figure  3-20  shows  the  peak  pointing  errors  that  serve  as  a 
function  of  the  amplitude  of  the  simulated  disturbance  torque  at  different 
RPM's.  (The  duration  of  the  disturbance  was  20  seconds.)  The  worst  case 
delays  errors  assumed  for  each  spin  speed,  and  the  available  motor  torque  and 
running  friction  were  assumed  to  be  170  in-oz  and  100  in-oz,  respectively. 

The  results  show  that  of  the  three  spin  speeds  tested,  the  system  is  least 
sensitive  to  disturbances  at  40  RPM.  It  should  also  be  noted  that  40  RPM  has 
the  best  limit  cycle  performance. 


^ ”*7t:n:  ^^9u*^e  3-10.  777  Search  Mode  with  Ground  Position  Loop  75  RPM 


= 4.8  sec 
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4.0  DESIGN  OF  GROUND  POSITION  AtjD  RATE  LOOP 

Subsequent  to  the  design  of  the  ground  position  loop,  it  was  desired  to 
study  the  effects  of  implementing  a rate  loop  in  the  ground  path  as  well,  by 
differencing  the  position  terms.  The  block  diagram  of  the  control  system 
with  both  a rate  term  and  a position  term  is  the  ground  loop  is  shown  in 
Figure  4-1 . 


4. 1 Linear  Analysis 

Unfortunately,  the  Z-plane  analysis  is  not  so  readily  done  in  this  case. 
That  is  because  the  rate  term  is  computed  by  differencing  position  measurements 
that  are  taken  at  a rate  of  once  per  second,  regardless  of  spin  speed.  However, 
the  spacecraft  receives  the  ground  command  only  once  per  sampling  period 
(Tg  = 2T)  as  before.  As  a first  attempt  to  analyze  the  problem, 
ihe  simplifying  assumption  will  be  made  that  the  rate  term  is  computed  by 
differencing  position  loop  measurements  taken  once  per  sampling  period.  Since 
this  is  actually  less  often,  i.e.,  worse,  than  what  is  actually  happening,  any 
results  obtained  by  this  analysis  will  be  conservative. 

The  open  loop  transfer  function,  including  delays, is  L(z)  where 


L(z)  = Lr(z)  + K^L^(z)  + K L2(z) 


(4-1) 


where 


L(^( z;  is  the  open  loop  transfer  function  of  the  rate  loop  in  the  spacecraft, 


Lr(z)  = K^'  KrK3 


2z(z-l) 


T 


(4-2) 


and  K^L^(2)  is  the  rate  term  in  the  ground  loop 


0 S X 


02  z + 

[ 2(z-l)^ 
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and  <1^(2)  is  the  position  term  in  the  ground  loop  or 


U2  z + z + Ujj 
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= .00017938  = 


^2  Kc*  Kn, 


= Selectable  gain 


= Ground  Rate  Gain 


Ground  Position  gain  = Kp'  x 261 
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Then  the  characteristic  equation  becomes 
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1 + L^(z)  + KpL2(z)  + *<^1-1(2^  = 0 


(4-5) 


In  order  to  determine  the  range  of  for  stability,  (K^  depends  on  the 

spin  period  T,  the  gain  state  chosen,  and  the  position  gain),  the  characteristic 

equation  will  be  rewritten  as 


KA,(t) 

1 + = 0 

1+4(2)  + KpL2(z) 


(4-6) 
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The  characteristic  equation  is 


1 > 


02  2-  + Oi  2+  Oq 

, P 0 s 

2 

W2  2 + Z + Uq 

22(z-1)^ 

X 

z 

6(z-l)^ 

^ 


02  Z + 


2(z-l)‘ 


= 0 


(4-7) 


3 k 

Multiplying  through  by  6(2-1)  2 , we  obtain 

zl'elz-l)^  + 3 z'"-''  [n2  z + 7 [z-l]  + KpK^*  [>^2  + w-,  z 


3 K 'K  K fz-1)  ion  ^ = 0 


(4-8) 


0 S X 


Putting  the  equation  into  the  form 


1 + G(z)  = 0. 


(4-9) 


we  have 
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Figures  4-2  through  4-4  plot  the  maximum  rate  gain  for  stability  as  a 
function  of  the  ground  loop  delay  for  various  gain  states.  The  graphs  are 
more  accurate  in  predicting  the  stability  of  the  actual  simulation  at  the 
higher  RPM's.  This  is  to  be  expected,  because  in  performing  the  analysis, 
the  assumption  was  made  that  the  rate  term  in  the  ground  loop  was  determined 
by  differencing  position  measurements  taken  once  per  sampling  period,  rather 
than  once  per  second,  as  is  done  in  the  simulation.  At  the  higher  RPM's, 
where  the  sampling  period  is  closer  to  1 second,  the  analysis  more  closely 
typifies  the  actual  flight  conditions,  whereas  at  the  lower  RPM's  (e.g.,  40) 
the  analytical  results  are  conservative.  In  any  case.  Figures  4-2  through 
4-4  are  to  be  used  only  as  a guide  to  determine  the  rate  gains,  and  the 
final  determination  must  be  made  by  simulation. 

4.2  Performance  Analysis 

The  performance  analysis  will  once  again  be  based  on  two  criteria  -- 
limit  cycle  performance,  and  the  transient  response  (peak  error)  due  to 
sudden  changes  in  the  friction  level.  The  results  will  show  that  there  is 
little  improvement  in  the  limit  cycles  behavior  when  the  rate  term  is 
included  in  the  ground  loop,  but  at  the  higher  RPM  there  is  some  improve- 
ment in  the  transient  response,  i.e.,  the  resulting  peak  errors  are  not 
as  large. 

Figure  4-5  shows  the  time  response  of  the  system  at  60  RPM  to  a simu- 
lated disturbance,  (10  in-oz  for  20  sec)  where  the  ground  loop  controller  has 
both  a rate  and  position  term.  As  can  be  seen,  the  limit  cycle  is  again  less 
than  2 degrees.  The  response  to  a larger  disturbance  (50  in-oz)  is  shown  in 
Figure  4-6.  Figure  4-7  summarizes  the  peak  error  responses  to  disturbance 
torques  of  various  amplitudes  at  60  RPM  and  compares  the  resulting  peaks  to 
the  maximum  error  when  the  ground  loop  controller  contained  only  a position 
term. 
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Figure  4-7. 
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At  75  RPM,  more  improvement  in  the  transient  response  is  shown  by  the 
inclusion  of  the  rate  term.  Typical  time  traces  are  shown  in  Figures  4-8 
and  4-9,  while  the  comparison  with  the  performance  of  the  controller  with 
a ground  position  loop  only  is  shown  in  Figure  4-10.  The  reason  that  more 
improvement  is  shown  at  the  higher  RPM  is  that  the  ground  rate  gain  is 
larger  compared  to  the  spacecraft  rate  loop  gain  at  75  RPM  than  at  60  RPM. 

At  40  RPM,  there  is  no  improvement  in  peak  performance  when  the  ground  i 

rate  loop  is  included,  and  that  is  because  the  ground  rate  loop  gain  is 
almost  insignificant  compared  to  the  rate  loop  gain  in  the  spacecraft. 

This  comparison  of  gains  is  summarized  in  Table  4-1.  ^ 

Table  4.1  Rate  Gains  ^ 


Figure  4-8.  Response  to  Small  Disturbance,  Ground  Position/Rate  Loop, 
75  RPM 
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^ Figure  4-10.  Comparison  of  Peak  Error  Response,  Ground  Position  Vs. 
Position/Rate  Loop,  75  RPM 
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5.0  IMPLEMENTATION 

The  implementation  of  the  ground  loop  controller,  i.e.,  the  criteria 
for  sending  a count  is  very  important,  because  the  quantization  neces- 
sarily produced  by  the  fact  that  only  an  integer  number  of  count  can 
be  sent  is  quite  large  compared  to  the  computed  control  command.  In  this 
report,  truncation  was  used  to  determine  when  commands  should  be  sent, 
while  round-off  is  to  be  used  in  the  actual  implementation  of  the  controller. 
This  is  the  reason  the  gains  were  divided  by  2 in  column  4 of  Table  5-1. 

In  addition,  1 count  in  either  direction  is  sufficient  for  controlling  the 
system,  so  the  counter  should  be  limited  to  ^ 1 count.  A diagram  of  the 
controller  implementation  is  shown  in  Figure  5-1. 

Some  question  has  been  raised  regarding  performance  when  the  plat- 
form has  been  knocked  off  the  earth,  and  the  telemetered  signals  are  no 
longer  valid,  but  instead  are  the  saturated  value.  Since  the  ground 
command  is  limited  to  1 count  anyway,  it  makes  no  difference  that  the 
values  telemetered  to  the  ground  controller  are  limited  to  + 6°  when  in 
reality  the  error  is  larger.  The  ground  controller  will  continue  to 
send  1 count  until  the  error  has  been  reduced  sufficiently  to  bring  the 
command  inside  the  deadzone. 

In  addition,  the  question  was  raised  as  to  how  long  the  error  should 
be  off  the  earth  before  the  control  system  is  inhibited.  In  the  normal 
mode,  the  criteria  is  3 counts,  or  3 seconds,  and  this  was  seen  to  be 
quite  insufficient.  It  is  recommended  that  the  criteria  of  60  counts 
or  1 minute  be  used  in  the  "fly  by  wire"  design,  in  order  to  allow  the 
control  system  more  than  sufficient  time  to  recover  from  large  transients. 
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TABLE  5-1  CONVERSION  FACTORS 
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6.0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  study  indicate  that  the  777  spacecraft  #9433 
can  be  controlled  with  a reasonable  degree  of  accuracy  by  operating 
the  spacecraft  in  the  search  mode  and  sending  up  commands  from  the 
ground.  It  should  be  realized  however,  that  because  of  the  large 
quantization  in  the  control  system,  a convergent  response  of  the 
pointing  error  can  never  be  obtained,  and  even  in  the  steady  state 
the  response  will  be  oscillatory. 

Two  basic  ground  loop  control  system  designs  were  considered. 

The  first  consisted  simply  of  proportional  control,  while  the  second 
used  proportional  plus  rate  control.  These  controller  designs  were 
evaluated  using  two  criteria  - limit  cycle  response  (steady  state) 
of  the  position  error,  and  the  transient  response  due  to  large  torque 
disturbances  (up  to  50  in  oz).  For  this  study,  the  designs  were  tested 
at  40,  60  and  75  RPM. 

In  the  steady  state,  the  amplitude  of  the  limit  cycle  oscillations 
depends  on  the  magnitude  of  the  position  gain,  and  the  stability  analysis 
shows  that  a higher  position  gain  is  possible  at  tl\e  lower  RPM's  (Figure 
3-1).  Using  the  gains  given  in  Table  5-1,  the  simulation  results  show, 
using  proportional  control  only,  that  both  the  best  limit  cycle  perform- 
ance and  transient  response  is  obtained  at  the  lower  RPM's.  At  40RPM,  a 
limit  cycle  of  about  1°  can  be  expected  compared  to  about  2 1/2°  at  75  RPM. 
The  peak  errors  in  response  to  transients  were  summarized  in  Figure  3-20. 
When  the  rate  loop  is  included,  some  improvement  is  seen  in  the  transient 
response  at  the  higher  RPM's  but  at  40  RPM  there  is  no  difference. 

It  is  thus  recommended  that  the  spacecraft  be  operated  at  the  lower 
RPM's.  In  addition,  it  is  recommended  that  the  control  system  first  be 
tested  using  only  the  position  loop,  and  if  this  is  satisfactory,  that 
the  rate  term  not  be  included. 
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APPENDIX  A.  GROUND  LOOP  TIME  DELAYS 

Reference:  DSCS-FX-004,  "Time  Delays  for  Ground  Loop  Control  System 
Development,"  W.  P.  Carey,  23  February  1976 

The  above  reference  contains  a detailed  analysis  of  the  various  time 
delays  associated  with  the  777  Comnand  and  Telemetry  lines  and  on-board 


despin  controller  measurement  logic. 

Table  A-1 , which 

is 

taken  from 

the  reference. 

presents  the  minimum  and  maximum  delays 

associated  with 

each  RPM,  The 

overall  minimum  value 

for  the  dairy  is 

2.4 

sec,  and  occurs 

at  75  RPM,  and 

the  maximum  is  6.4  sec 

, occurring  between  ■ 

47  and  50  RPM. 

The  delays 

are  plotted  as  a function  of  spin  rate 

in 

Figure  A-1 . 

Table  A-1 

Time 
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RPM 

T 
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Tp  fiAX 

75 

.8 

2.4 

4.8 

70.4 

.852 

Not  Given 

5.11 

70 

.857 

2.57 

14.29 

65 

.923 

2.77 

4.62 

60 

1.0 

3.0 

5.0 

55 

1 .091 

3.27 

5.46 

50 

1.2 

3.6 

6.0 

46.97 

1.279 

Not  Given 

5.11 

45 

1.333 

4.0 

5.33 

44.7 

1.342 

2.64 

Not  Given 
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APPENDIX  B.  SIMULATION  DESCRIPTION 

The  simulation  results  presented  in  this  document  were  obtained  from 
the  Tymshare  dig-ital  simulation  briefly  described  in  this  appendix.  This 
simulation  (as  well  as  this  description)  is  based  on  the  same  111  simulation 
used  in  Reference  3,  with  some  modifications.  The  capabilities  and  limit- 
ations of  the  simulation  are  discussed  in  the  following  paragraphs,  as  well 
as  the  main  differences  between  this  simulation  and  the  one  of  Reference  3. 

B-1  Simulation  Model,  Saturation  Limits 

The  simulation  model  is  shown  in  Figure  B-1.  The  simulation  is  patterned 
closely  after  the  hardware  limitations  in  order  to  realistically  evaluate  the 
effect  of  hardware  non-linearities.  For  example,  each  of  the  saturations 
shown  in  the  simulation  represent  saturation  levels  in  the  hardware.  These 
saturation  limits  were  often  reached  when  the  simulation  was  operated  in 
the  search  mode.  However,  in  the  "fly  by  wire"  mode  of  operation,  where 
at  most  a few  counts  are  being  cownanded  from  the  ground  loop  and  space- 
craft rate  loop  combined,  these  can  be  neglected  for  the  most  part.  For 
example  the  number  counts  needed  to  reach  some  of  the  saturation  limits  is 
si'iown  below: 

I-j  ( = 14  volts) 

k2k^Hi>  14  or  H^>  1942  counts,  for  saturation 
I2  ( = 14  volts) 

*^1*^9^1^  14  or  H.j>  14  counts,  for  saturation 

A more  complicated  model  of  the  integrator  in  the  rate  mode  (denoted  by 
levels  Ig,  Ig)  showing  the  saturation  effects  was  introduced  in  Reference  4, 
and  this  model  has  been  retained  in  the  simulation,  although  only  the  behavior 
in  the  linear  region  is  indicated  in  the  diagram. 

B-2  Rate  Measurement  Logic  and  Relay 

The  measurement  logic  is  a functional  representation  of  the  hardware. 

The  platform  rate  error  signals  are  not  obtained  in  this  simulation  as  they 
are  in  the  hardware,  by  processing  the  earth  sensor  and  platform  reference 
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pipper  signals.  Instead,  the  up-down  counter  contents  are  calculated  from 
the  platform's  instantaneous  rate  using  appropriate  scale  factors. 

The  model  of  the  measurement  logic  in  the  rate  loop  of  the  spacecraft 
is  given  in  Figure  B-2.  This  model  is  different  than  that  of  the  simu- 
lation of  Reference  3,  which  was  incorrectly  modeled  witli  a double  sided 
deadband,  as  in  Figure  B-3a.  (For  the  mode  of  operation  for  which  Refer- 
ence 3's  simulation  was  intended,  this  inaccuracy  is  insignificant.)  For 
the  "fly-by-wire"  design,  when  the  rate  loop  is  at  most  putting  out  a few 
counts,  and  very  often  only  one  count,  the  modeling  of  the  counts  around 
zero  is  very  significant.  A provision  has  been  made  to  study  the  effects 
of  an  asymmetric  quantization  model,  (see  Figure  B-2)  since  in  reality  the 
counter  is  not  "dead  set"  around  0. 

The  measurement  delay  of  one  quarter  of  a sampling  period  in  the  space- 
craft rate-loop  is  implemented  by  saving  the  necessary  variables  3/4  of  the 
way  through  the  integration  subroutine,  to  be  used  in  the  next  sampling 
period. 

B-3  Ground  Loop  Model 

The  gound  loop  consists  of  both  a position  and  rate  term.  The  position 
gain  Al.  is  scaled  the  same  way  as  Kp*"  in  the  analysis  presented  in  the 
body  of  the  report.  That  is,  the  actual  position  gain  will  be  4.56  Al 
counts/degree.  This  scaling  can  be  changed  by  amending  line  845  if  the  user 
finds  some  other  units  more  convenient. 

The  rate  term  is  computed  by  differencing  position  terms  taken  once 
per  second,  regardless  of  spin  speed.  The  rate  gain,  A4,  is  in  units  of 
counts/radian/sec,  but  this  can  also  be  changed  by  altering  line  843. 


Counts 


Rate  (Counts) 


Figure  B-3a.  Incorrect  Model  - 
Double  Deadband  at  Origin 


Counts 
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The  method  of  quantization  (truncation)  has  already  been  indicated  in 
Section  2.0,  and  one  count  is  sent  when  the  magnitude  of  the  ground  command 
exceeds  1.  If  it  is  desired  to  change  this  to  round-off,  the  position  and 
rate  gains  presented  in  the  analysis  should  be  halved. 

B-4.  Implementation  of  Ground  Loop  Delay 

When  this  study  of  the  “fly-by-wire"  design  first  began,  and  only  a 
position  loop  was  considered,  the  ground  loop  delay  was  simply  modeled  as 
an  integer  number  of  sampling  periods,  or  = 06  T^,  where  06  is  an  input 

to  the  program.  (If  06  = 0,  the  ground  loop  delay  was  that  of  the  rate  loop, 
or  T^/4.)  The  following  diagram  shows  how  this  was  done. 


2T. 


3T 


+ 


3T 


y^ 


X 


V 


Figure  B-4.  Implementation  of  Ground  Loop  Relay  in  Simulation 
(Position  Loop) 

The  fact  that  XI  was  sampled  "immediately  before"  T^  is  indicated  by  the  use 
of  the  variable  XI  in  line  845.  The  range  of  values  of  06  to  be  used  were 
taken  from  Appendix  A.  The  fact  that  position  measurements  were  taken  only 
once  per  second  was  not  modeled,  since  any  delays  associated  with  this  could 
be  accounted  for  by  an  appropriate  choice  of  06. 

When  the  rate  term  in  the  ground  loop  was  added,  it  was  necessary  to 
model  the  fact  that  the  samples  were  taken  only  once  per  second  since  the 
rate  term  is  computed  by  differencing  position  measurements.  Figure  B-5 
shows  the  implementation  of  the  delay  under  those  conditions. 
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1 sec  2 sec  3 sec  etc. 


Figure  B-5.  Implementation  of  Ground  Loop  Delay  - Measurements  Taken 
Once/Sec 

The  important  point  to  notice  is  that  the  resulting  delay  will  be  greater 
than  D6  T^.  (For  many  of  the  simulation  results  the  position  term  in  the 
ground  loop  was  still  taken  as  shown  in  Figure  B-4,  i.e.,  for  the  position 
term  only,  the  fact  that  measurements  were  taken  only  once  per  second  was  not 
used.  This  can  be  easily  fixed  by  replacing  XI  by  X(l)  in  line  845,  and 
any  future  simulation  work  should  be  done  with  this  change  made.) 

B-5.  Dynamics 

The  platform  dynamics  are  simply  modelled  in  the  Tymshare  simulation  as 
rigid  body,  single  axis  linear  dynamics,  i.e.  there  is  no  cross  coupling  with 
the  rotor  or  other  platform  rates.  Nutation  cross-coupling  and  damper 
dynamics  are  ignored.  However,  comparison  between  the  dynamics  of  this  simu- 
lation at  the  time  it  was  originally  written,  and  design  verification  runs, 
indicate  that  satisfactory  results  were  obtained  from  using  the  simplified 
model . 

B-6.  Motor  Model 

Dne  aspect  of  the  simulation  model  that  has  recently  come  under  close 
scrutiny  is  the  motor  model.  Presently,  the  maximum  motor  capability  is 
modeled  as  being  constant,  or  the  same  for  all  spin  speeds.  However, 
recent  studies  have  indicated  that  the  motor  capability  might  be  greatly 
reduced  at  the  higher  RPM's.  At  this  time  no  more  definite  model  has  been 
presented,  and  so  this  limitation  has  not  been  taken  into  account  in  the 
simulation. 


TABLE  B-1 

PARAMETER  DEFINITIONS  FOR  PROGRAM  777 's 
"FLY-BY-WIRE"  TYMSHARE  SIMULATION 
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PROGRAM  NOMINAL 

SYMBOL  value 


DEFINITION 


F 

1 

PI 

1 

N7 

20 

XI 

variable 

X2 

variable 

X4 

variable 

X5 

360 

X6 

variable 

B(l) 

0 

B(3) 

0 

R 

0 

NO 

0 

N1 

variable 

11 

14 

12 

14 

15 

5 

16 

-14 

B3 

.0667 

FI 

1643 

K1 

5.86x10'^ 

K2 

1.23 

K3 

.131 

K4 

.2286 

K5 

4.48 

K6 

variable 

Mode  control 

F » 1 - Inhibit  transfer  from  rate/search 
mode  to  normal  control 

Number  of  sample  periods  between  printout 

Number  of  integration  steps  per  sampling  interval 

Platform  position  error  (deg) 

Platform  rate  error  (deg/sec) 

Integrator  voltage  in  lead-lag  filter  (volts) 

Spacecraft  rotor  spin  rate  (deg/ sec) 

Relative  rate  between  platform  and  rotor 
(deg/sec)  X6  = X5  - X2 

Measurement  voltage  offset  during  rate/?earch 
mode 

Platform  rate  bias  during  serach  mode  (counts) 

Spacecraft  rate  error  (deg) 

Initialization  of  random  number  sequence 

Measurement  noise  during  serach  mode 

Saturation  level  (volts) 

Saturation  level  (volts) 

Saturation  level  (volts) 

Saturation  level  (volts) 

Integrator  gain  during  rate/search  mode 

Measurement  clock  frequency  during  rate/ search 
mode  (counts/sec) 

D/A  converter  gain  (volts/count) 

Gain  during  rate/search  mode  which  replaces 
the  lead/lag  network  (volts/volt) 

Integrator  gain  (volt/volt) 

Relative  integrator  gain  (volt/volt) 

Relative  lead-lag  gain  (volt/volt) 

Ground  commandable  controller  gain  (volt/volt) 


\ 

\ 


■i 
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PROGRAM 

SYMBOL 

NOMINAL 

VALUE 

DEFINITION 

K7 

.4 

Motor  gain  (ft-lb/volt) 

K8 

1 

Measurement  scaling  gain  (not  really  used  In 
simulation) 

K9 

183.9 

Integrator  Input  path  gain  (volt/volt) 

MO 

170 

Motor  saturation  (In-oz) 

Ml 

100 

Motor  friction  (In-oz) 

M2 

.14  X6 

Motor  back  emf  (In-oz) 

J 

69 

Platform  pitch  Inertia  (s1ug-ft^) 

J1 

339 

Rotor  pitch  Inertia  (s1ug-ft^) 

06 

variable 

Delay  In  ground  loop  (number  of  sampling  periods) 
(see  Appendix  A)  (If  D6  = D,  Delay  Is  same  as 
spacecraft  rate  loop  delay,  I.e.,  .5T) 

A1 

.15 

Ground  loop  position  gain  (K') 

A4 

130 

Ground  loop  rate  gain,  counts/radlan/sec 

LI 

1 

Ground  loop  counter  limit 

A5 

.5 

Offset  bias  In  simulation  spacecraft  rate  loop 
counter  (see  Figure  B-2) 

TO 

variable 

Simulation  run  time  (sec) 

s(i) 

variable 

Sample  period  during  rate/search  mode 

M9 

variable 

Number  of  torque  pulses  to  be  simulated 

D(I) 

variable 

Amplitude  of  1-th  pulse  (In-oz) 

U(I) 

variable 

Time  1-th  pulse  begins  (sec) 

G(I) 

variable 

Time  1-th  pulse  ends  (sec) 

H 


1 
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